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We present a theoretical discussion of a self-assembly scheme which makes it possible to use 
DNA to uniquely encode the composition and structure of micro- and nanoparticle clusters. These 
anisotropic DNA-decorated clusters can be further used as building blocks for hierarchical self- 
assembly of larger structures. We address several important aspects of possible experimental 
implementation of the proposed scheme: the competition between different types of clusters in a 
solution, possible jamming in an unwanted configuration, and the degeneracy due to symmetry with 
respect to particle permutations. 
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Over the past decade, a number of proposals have iden- 
tified potential applications of DNA for self-assembly of 
micro- and nonostructures [H,3,[35,3,[55. Among these 
proposals, one common theme is finding a way to utilize 
the high degree of selectivity present in DNA-mediated 
interactions. An exciting and potentially promising ap- 
plication of these ideas is to use DNA-mediated interac- 
tions to programmable self-assemble nanoparticle struc- 
tures , [7] , [8] , [9] . Generically, these schemes utilize 
colloidal particles functionalized with specially designed 
ssDNA (markers), whose sequence defines the particle 
type. Selective, type-dependent interactions can then 
be introduced either by making the markers complemen- 
tary to each other, or by using linker-DNA chains whose 
ends are complementary to particular maker sequences. 
Independent of these studies, there are numerous pro- 
posals to make sophisticated nano-blocks which can be 
used for hierarchical self-assembly. One recent advance 
in the self-assembly of anisotropic clusters is the work of 
Manoharan et. al[10]. They devised a scheme to pro- 
duce stable clusters of n polystyrene microspheres. The 
clusters were assembled in a colloidal system consisting 
of evaporating oil droplets suspended in water, with the 
microspheres attached to the droplet interface. The re- 
sulting clusters, unique for each n, are optimal in the 
sense that they minimize the second moment of the mass 
distribution M2 = YJi^i^^i - ^cmY- 

In this paper, we present a theoretical discussion of 
a method which essentially merges the two approaches. 
We propose to utilize DNA to self-assemble colloidal clus- 
ters, somewhat similar to those in Ref. [10]. An im- 
portant new aspect of the scheme is that the clusters 
are "decorated": each particle in the resulting cluster is 
distinguished by a unique DNA marker sequence. As a 
result, the clusters have additional degrees of freedom as- 
sociated with particle permutation, and potentially may 
have more selective and sophisticated inter-cluster inter- 
actions essential for hierarchic self-assembly. In addi- 
tion, the formation of such clusters would be an impor- 
tant step towards programmable self-assembly of micro- 
and nanostructures of an arbitrary shape, as suggested 
in Ref. [TH. 

We begin with octopus-like particles functionalized 
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FIG. 1: (Golor online) The minimal second moment clusters 
for n — 5, 6, 7, and 9. Pictures of all the clusters from n — A 
to 15 are available in p!o| . 



with dsDNA, with each strand terminated by a short 
ssDNA marker sequence. We assume that each parti- 
cle i has a unique code, i.e. the maker sequence Si of 
ssDNA attached to it. We then introduce anchor DNA 
to the system, ssDNA with sequence s^SB-.-Sn, with Si 
the sequence complementary to the marker sequence Si. 
The anchor is designed to hybridize with one particle of 
each type. Consider a cluster of n particles attached to 
a single anchor. 

If we treat the DNA which link the particles to the 
anchor as Gaussian chains, there is an entropic contribu- 
tion to the cluster free energy which can be expressed in 
terms of the particle configuration {ri, ...,r^} as follows. 
Here Rg is the radius of gyration of the octopus-like DNA 
arms. 
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This approximation of the DNA arms as Gaussian chains 
is acceptable provided their length L exceeds the per- 
sistence length Ip 50 nm and the probability of self- 
crossing is small[12]. The physical mechanism which 
determines the final particle configuration in our system 
is quite different from the capillary forces of Manoharan 
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FIG. 2: (Color online) Schematic representation of the 
method for constructing decorated colloidal clusters using ss- 
DNA "anchors". 



et al. However, because the functional form of the free 
energy is equivalent to the second moment of the mass 
distribution, the ground state of the cluster should cor- 
respond to the same optimal configuration. 

Consider a system with n particle species and an an- 
chor of type saSB'-'Su' The clusters we would like to 
build contain n distinct particles (each particle in the 
cluster carries a different DNA marker sequence) at- 
tached to a single anchor. Let Cn denote the molar 
concentration of the desired one anchor cluster. Because 
there are many DNA attached to each particle, multiple 
anchor structures can also form. The question is whether 
the experiment can be performed in a regime where the 
desired one anchor structure dominates, avoiding gela- 
tion. 

We consider the stability of type Cn with respect to 
alternative two anchor structures. To do so we de- 
termine the concentration Cn+i of n + 1 particle struc- 
tures which are maximally connected, but do not have a 
1 : 1 : • • • : 1 composition. In particular, these structures 
contain more than one particle of each type, which could 
cause problems in our self-assembly scheme There 
are also n particle structures Cn with the correct compo- 
sition, but which contain two anchors. We would like to 
avoid the formation of these structures as well, as their 
presence decreases the overall yield of type Cn- Figure 
[3] enumerates the various structures for an n = 3 species 
system. If the experiment can be performed as hoped, 
we will find a regime where the ratios and ^ are 

small. To this end, the equilibrium concentrations Cn, 
Cn+i, and Cn are determined by equilibrating the chem- 
ical potential of the clusters with their constituents. 

Let Ci denote the molar concentration of species i, and 
Ca the molar anchor concentration. We consider the 
symmetrical case AGi = AG and equal initial particle 



concentrations cf"^ = c^^^ for all species i. 
we have ca = cb = = Cn = c. 



In this case 
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FIG. 3: (Color online) The topologically distinct one and 
two anchor structures for an anchor ssDNA with sequence 
s AS BSC- Different structure varieites may be obtained by 
relabeling the particle indices subject to the constraint that 
no more than one particle of each type is attached to a given 
anchor. 
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The binding free energy of the cluster has a contri- 
bution from the hybridization free energy AG associ- 
ated with attaching a particle to the anchor, and an en- 
tropic contribution from the number of ways to construct 
the cluster (since each particle has N hybridizable DNA 
arms). In addition we must take into account the en- 
tropy for the internal degrees of freedom in the structure 
stemming from the flexibility of the DNA attachments to 
the anchor. In the Gaussian approximation, neglecting 
the excluded volume between particles this localization 

volume V = (^) ^ i?^ can be calculated exactly by inte- 
grating over the various particle configurations weighted 

by the Boltzmann factor: = J (firi...(Prn exp 

We now consider the competing two anchor struc- 
tures Gn+i and Cn- The localization volumes V2 = 

-3 -3(n-l) -3 -3n 

2(-+2) V and ^3 = n^c^+i) 2 2(-+i) ^ can be calcu- 
lated in a similar fashion to v. Since there are many 
DNA attached to each particle, in what follows we omit 
factors of ^^^7^. 
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The concentration of free anchors Ca can be determined 
from the equation for anchor conservation. 



: Ca + Gn + 2Gn + 2nGn+l 



(5) 
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FIG. 4: (Color online) The molar concentrations Ca, Cn, 
nCn+i, and Cn in the symmetrical case for a system with 
n = 5 particle species. The total particle volume fraction 
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We are interested in the low temperature regime 
where there are no free anchors in solution. We de- 
termine the saturation values for the ratios of interest by 



noting that the Boltzmann factor S = exp 
in this regime. 
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Since 



- <C 1, eq. [71 provides the experimental con- 
straint for suppressing the two anchor structures. Taking 
the radius of the hard spheres R ^ Rg^ it can be inter- 
preted as a criterion for choosing the initial anchor con- 
centration ci^^ for an n species system with (j) = ^Rg 
the particle volume fraction for an individual species. 
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The condition gives the maximum anchor concentration 
for the two anchor structures to be suppressed. Since 



< - the theoretical limits are 



< 1, .29, .06, and .01 



for n = 4, 5, 6, and 7 respectively. In figure [H we plot the 
solution for the concentrations. There is a large temper- 
ature regime(^^ < 2) where the two anchor structures 
are suppressed in favor of the desired one anchor struc- 
tures. 

The experimental criterion for suppression of the two 
anchor structures (eq. [8|) provides a fairly strict bound on 
the anchor concentration, and hence the cluster yield for 
n > 5. However, the previous discussion considers only 



the equilibrium concentrations. Below the melting tem- 
perature (see Figured)), the connections are nearly irre- 
versible. In the irreversible regime where the DNA bind- 
ing is very strong, the probability that a cluster of type 

Cn becomes attached to a second anchor is p ~ ^^oj- 
The factor of n counts the possible attachment sites to 

the cluster. In terms of the fraction fn ~ of parti- 
cles in clusters of type Cn, the fraction of particles /2 in 

the two anchor structures will be /2 ~ |/n = 

Here the factor of \ is necessary since the initial attach- 
ment can be to either of the two anchors. As indicated 
in Figure [5l provided the anchor concentration is dilute 
enough, > 4, the formation of the two anchor struc- 
tures is suppressed. 
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FIG. 5: (Color online) Cluster fractions as a function of the 

concentration ratio for an n = 5 species system in the 

irreversible binding regime, fn is the fraction of the desired 
one anchor structures and /2 is the fraction of two anchor 
structures. 

We now present a brief discussion of the role that jam- 
ming plays in preventing the one anchor structures from 
assuming the minimal second moment configuration. We 
performed simulations of the assembly of optimal col- 
loidal clusters up to n = 9 particles by numerically inte- 
grating the particles' Langevin equations. 



dt 



-ViH + rii 



(9) 



Here b is the particle mobility, and the thermal noise 
has been artificially suppressed(i.e. 77 = 0). The model 
Hamiltonian H used in these simulations has been dis- 
cussed in detail elsewhere pj]. As indicated in Figure 
[6l the hard sphere system gets trapped in a configura- 
tion with a larger M2 than the optimal cluster, whereas 
the soft-core system is able to fully relax. The jam- 
ming behavior is largely determined by the single control 
parameter with d the diameter of the hard sphere. 

Beyond the critical value ^ ^ -5 the jamming behavior 
is either completely eliminated, or greatly reduced in the 
case of larger clusters. 
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FIG. 6: (Color online) Plot of the dimensionless second mo- 
ment ^ as function of time for n — 9 particles. Results are 
shown for the case of hard spheres and also for a system with 
a soft-core repulsion with geometric parameter = .7. The 
dashed line is the theoretical moment for the triaugmented 
triangular prism, which is the minimal n — 9 structure. 

Building these decorated colloidal clusters is the first 
major experimental step in a new self-assembly proposal. 
However, in order to utilize the resulting clusters as build- 
ing blocks, an additional ordering is necessary. The 
problem is that the decoration introduces degeneracy in 
the ground state configuration. This degeneracy was 
not present in [10] since all the polystyrene spheres were 
identical. Namely, in the colloidal clusters self-assembled 
by our method, permuting the particle labels in a cluster 
does not change the second moment of the mass distribu- 
tion (see Figure [7j). We need a method to select a single 
"isomer" out of the many present after self-assembly. In 
the DNA-colloidal system considered here, this isomer se- 
lection can be facilitated by "linker" ssDNA. These are 
short ssDNA with sequence sasb to connect particles A 
and B. We first construct a list of nearest neighbors for 
the chosen isomer, and introduce linker DNA for each 
nearest neighbor pair. The octopus-like DNA arms of 
the given particles will hybridize to the linkers, resulting 
in a sping-like attraction between the selected particle 
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pairs. Note that the length L of the DNA arms must be 
on the order of the linear dimension of the original clus- 
ter. Otherwise the interparticle links cannot form upon 
introduction of linker DNA to the system. It should be 
noted that although this method breaks the permutation 
degeneracy of a cluster, the right-left degeneracy will still 
be present. 

In conclusion we discussed a method which uses DNA 
to self-assemble anisotropic colloidal building blocks. We 
found an experimentally accessible regime where the re- 
sulting clusters are minimal second moment configura- 
tions. In addition, the clusters are decorated: each par- 
ticle in the cluster is distinguished by a unique marker 
DNA sequence. The cluster formation process pro- 
vides an interesting model system to study a new type of 
jamming-unjamming transitions in colloids. Construct- 




FIG. 7: (Color online) An illustration of degeneracy in DNA- 
coded nanoclusters. Two different n = 6 isomers are pic- 
tured, both with the same minimal second moment configu- 
ration, the octahedron. 



ing decorated colloidal clusters would represent a major 
step towards realizing the long-term potential of DNA- 
based self-assembly schemes. 

Acknowledgement. This work was supported by 
the ACS Petroleum Research Fund (grant PRF #44181- 
ACIO), and by the Michigan Center for Theoretical 
Physics (award # MCTP 06-03). 



[8] Winfree, E.; Liu, F.; Wenzler, L. A.; Seeman, N. C. 

Nature 1998, 394, 539-544. 
[9] Cobbe, S.; Connolly, S.; Ryan, D.; Nagle, L.; Eritja, R.; 

Fitzmaurice, D. J. Phys. Chem. B 2003, 107, 470-477. 
[10] Manoharan, V. N.; Elsesser, M. T.; Pine, D. J. Science 

2003, 301, 483-487. 
[11] Licata, N. A.; Tkachenko, A. V. arXiv.cond- 

mat/05044 07 v2 2006, . 
[12] Tkachenko, A. V. Physical Review Letters 2002, 89, 

148303. 



